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I. INTRODUCTION 

Fresnel diffraction at a straight edge leads to an inten- 
sity pattern determined by the Cornu spiral [l|, The 
latter results from an integral with a quadratic phase. 
The Talbot effect Q is a generalization of this scattering 
situation. The discreteness of the grating translates the 
integral of quadratic phase factors into a sum. Recently, 
the combination of chirped pulses together with an ap- 
propriate atomic level scheme created an atomic analog 
of the Fresnel diffraction of the straight edge0, 01 • In 
the present paper we propose the atomic analog of the 
Talbot effect by chirping a two-photon transition in a 
multi-state ladder. 

Chirped pulses are characterized by a non-linear phase 
dependence and have a vast variety of applications [1, 0, 

EMM 

.[T]|. Examples include, the realization of a time- 
domain Fresnel lens with coherent control andquan- 
tum state measurement using coherent transients [ISj- Of 
particular relevance in the context of the present pa- 
per are interference phenomena in quantum ladder sys- 
tems induced by two-photon excitation with chirped laser 
pulses[il[i|]. 

The central idea is summarized in Fig. [TJ We consider 
a two-photon transition from the ground state \g) to the 
excited state |e). In the case of a single intermediate 
state[l^ discussed in the top row of this figure the 
total excitation probability We = |cep is due to the inter- 
ference of two paths: (?) in the direct path two photons 
of energy hioo are absorbed instantaneously and (ii) in 
the sequential excitation path the total energy of 2huJo 
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is absorbed in two steps: in the first step the system is 
excited to the intermediate state |to), whereas the sec- 
ond, delayed step provides the residual energy to reach 
the excited state |e). This interference effect gives rise 
to an excitation probability determined by a Fresnel in- 
tegral. It is experimentally accessible by measuring the 
population in the excited state through the detection of 
the fluorescence signal. Population transfer with chirped 
laser pulses has been experimentally demonstrated in the 
three-state ladder of Rubidium Q . 

In the case of a multi-state ladder with a manifold of 
intermediate states shown in the second row of Fig. [T] we 
find a generalization of this two-path interference to a 
multi-path interference. Here the interference is mainly 
due to the indirect paths through the individual interme- 
diate states. The direct path does not play a dominant 
role anymore. The problem of scattering a chirped pulse 
from D intermediate levels is very much in the spirit of 
diffraction from a Z?-slit grating. For an equidistant man- 
ifold the phase of the excitation path through the m-th 
level contributing to the total excitation probability de- 
pends quadratically on m in complete analogy to the 
Talbot effect. We emphasize that the present proposal 
is different from the temporal Talbot effect suggested in 
Ref. [11. 

Our paper is organized as follows: In Sec. |IT1 we in- 
vestigate the population transfer in the D-state ladder 
system shown in Fig. [1] driven by a weak chirped laser 
pulse. We use second order perturbation theory to de- 
rive an analytic expression for the probability amplitude 
to be in the excited state. This result is of the form 
of a Gauss sumfiol. [20| which plays a prominent role in 
number theory. We dedicate the following sections to a 
discussion of the physical origin of this Gauss sum. 

We start in Sec. IIIII bv first considering the transition 
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probability amplitude Ce in the limit of a single inter- 
mediate state. In this case Ce follows from the comple- 
mentary error function with a complex-valued argument 
which depends on the chirp parameter. With the help 
of appropriate asymptotic expansions of the complemen- 
tary error function we bring to light the oscillations in 
the transition probability We = |cep shown in Fig. [T] for 
negative values of the chirp. These oscillations arise due 
to the interference of two excitation paths as shown in 
Sec. I IV I Moreover, a primitive version of the method of 
stationary phase allows us to identify the times of the 
transitions. We show that the accumulated phases de- 
pend quadratically on the offset Sm- Since our analy- 
sis works in frequency-time pha se space our approach is 
different from Refs. [3, ITsl. l2ll| which relies on the fre- 
quency domain. In the case of many intermediate levels 
the interference of quadratic phase factors gives rise to 
a Gauss sum which we express in its canonical form in 
Sec. El Here all dimensionless quantities are related to 
experimental parameters. In Sec. IVII we conclude by ad- 
dressing possible experimental realizations of multi-state 
systems. 



II. EXCITATION PROBABILITY 

We consider the multi-state ladder system shown in 
the bottom of Fig.[T]driven by a weak chirped laser pulse 
and calculate the probability amplitude of the excited 
state in second order perturbation theory. We show that 
this probability amplitude is the interference of many 
amplitudes which depend quadratically on the offset of 
the intermediate state. We present a simple phase space 
argument for the origin of such quadratic phase factors. 



A. Model system 

The ladder system consists of a ground state \g), an ex- 
cited state |e) separated by an energy 2hu;o and interme- 
diate states |m) with quantum numbers M' < m < M , 
as depicted in Fig. [1] These intermediate states are dis- 
placed by the offset 



6 + m A. 



(1) 



with respect to the central frequency loq. The harmonic 
manifold is characterized by the offset S oi a, specific cen- 
tral state |0) and separation A of adjacent states. 

In the interaction picture the Hamiltonian describing 
the interaction between the electric field of the chirped 
pulse 
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FIG. 1: Chirping a two-photon excitation through a single 
intermediate state (top) and through an equidistant manifold 
(bottom). The left and right columns show the level system 
and the resulting excitation probabilities as a function of the 
chirp parameter <f)" , respectively. In both systems the ground 
state \g) is connected by a two-photon transition to the ex- 
cited state je). In the case of a single intermediate state |m) 
with the offset S > with respect to the central frequency 



LJo discussed in the top row the population 



of the ex- 



cited state displays oscillations for negative chirps, whereas 
for positive values of the chirp it is decreasing rapidly. The 
oscillations result from an interference of the direct and the 
indirect excitation. The direct path is a direct two-photon 
transition through the virtual state of frequency loq. The in- 
direct path consists of first an excitation to the level |m) and 
then a transition to |e). Here we have used the parameters 
5 = 0.0225 fs"S Aio = 0.1525 fs"\ Qem^mg = 1 as suggested 
by the experiment . In the bottom row we discuss the 
case of a manifold of M' + AI + 1 equidistant intermediate 
states |m) with — Af < m < M which are shifted by the 
offset Sm = S + mA with respect to the central frequency 
ujo- The offset of the central state is S whereas the states in 
the harmonic manifold are separated by A. The population 
|cep of the excited state is due to interference of compet- 
ing quantum paths via the intermediate states and displays 
strong modulations as a function of the chirp. Parameters 
are M' = 7, M = 7, S = 0.0225 fs'S A = 0.003 fs'S Au = 

0.1525 fs , i^em^mq 1 



function f{t) and our ladder system reads 
V{t) = -ft^ [ans e'^'"*/(i) \m){g\ + 



ne,ne-'^"^'f{t)\e){m\+c.c. 



(3) 



of amplitude £q, carrier frequency loq and pulse shape 



where we have introduced the Rabi frequencies Qjj = 
Sopij/h with the dipole matrix element pj„i of the tran- 
sition \j) |m) and \ j) = \g,e). 



3 



B. Time evolution 



C. Description of the chirped pulse 



Throughout the paper we consider a situation in which 
the laser pulse is short (femtosecond-time scale) com- 
pared to the characteristic time scale of radiative decay 
processes. As a consequence we can neglect spontaneous 
emission and the time-evolution of the state vector 



k 

is given by the Schrodinger equation 



(4) 



(5) 



Since the laser field is assumed to be weak, the prob- 
ability amplitude Cg of the excited state can be simply 
calculated in second order perturbation theory. We find 



Ce{t) 



1 

2ff 



t t' 

dt' j dt" {e\V{t')V{t")\g) (6) 

oo — oo 



The first-order term vanishes since we assume that at 
to = — oo, that is before the interaction, the ladder sys- 
tem is initially prepared in the ground state \g) and a 
two-photon transition is necessary in order to reach the 
excited level. 

When we substitute the explicit expressions for the in- 
teraction Hamiltonian, Eq. ([3]), and concentrate on times 
after the pulse has passed, that is t — > oo, we obtain 



(7) 



for the probability amplitude Cg = Ce(t — > oo). 
The domain of the integral 

oo t 

X,„= / dt' j dt"e-* *-"(*'-*")/(*')/(*") (8) 



-oo — oo 



extends over the area below the diagonal crossing the 
(t', t")-plane from the lower left to the upper right. When 
we introduce the new integration variables i=t' + 1" and 
t = t' — t" and use the formula 



oo oo 



dt' / dt" • • • = - / dt dt. 



(9) 



-oo — oo 



-OO 



we eventually arrive at 



4 

m 

oo oo 

di j dt e~ 

oo 



t + t 



t-t 



(10) 



This expression is valid for an arbitrary pulse shape /. In 
the next section we restrict ourselves to a linearly chirped 
pulse with a Gaussian envelope. 



We now specify the pulse shape 



f(t) = fo cxp 
with the complex-valued amplitude 



^{^cofoft^ 



1 + ia 



(11) 



(12) 



and the dimcnsionlcss parameter a denotes the second 
order dispersion 



2 J," 



a = Auj (j)' 



(13) 



Here Alj is the bandwidth of the pulse and ((>" ~ 
d^ (l>{uj) I doj"^ is a measure for the quadratic frequency de- 
pendence of the phase of the laser pulse. 
It is useful to decompose the argument 



^(Atj/o)^ = Ur + iai 



(14) 



of the Gaussian pulse shape given by Eq. (fTTj) into the 
real and imaginary parts 



and ai 



(15) 



2 1 + 2 H-o2 

resulting in 

f{t) ^ fo cxp {-art^) cxp (-i a^t^) . (16) 

This representation brings out most clearly that the elec- 
tric field Eq. of the chirped laser pulse features a linear 
variation of the instantaneous frequency 



U!(t) ^ ujQ + 2ai t. 



D. Excitation probability amplitude 



(17) 



When we now substitute the pulse shape, Eq. (|11|) . 
into the expression Eq. (fTUl) for the excitation probability 
amplitude we obtain 



with 



oo 

h^fo f dte-ii^-for^ = '-^ 
J Auj 



(18) 



(19) 



and 



h = k j dt e-'«'"*e-3(A"/o) t . (20) 
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The substitution 



Aw 5„ 



2 ^ Aujfo 

casts the integral I2 into the form 



Vm 

fo 



r I ■ Vm 

eric I—— 



where we have introduced the dimensionless offset 

Vm — » 

and have recalled[22l the definition 



erfc(C) = —;= dz e 

V7I" J 

C(C) 



(21) 



(22) 



(23) 



(24) 



of the complementary error function. 

The integration in Eq. (PH)) over t from to cx) trans- 
lates into a path C(Cm) in the complex plane which ac- 
cording to Eq. ([21]l starts at the point Cm = i Vm/ fo and 
follows a straight line to infinity. This path encloses an 
angle arg(/o) with respect to the positive real axis. 

When we introduce the abbreviation 



the definition Eq. of fo yields the expression 
Ce = ^ w„i exp (i rjl^a) 



(25) 



(26) 



for the transition probability amplitude Eq. (|T8)l . Here 
we have introduced the weight factors 



-^erfc(C™(a))exp (-77^ 



with 



Cm (a) = Vm i Vl - «a. 



(27) 



(28) 



Equation (p6)) is the central equation of the present pa- 
per. It demonstrates that the probability amplitude Ce 
to be in the excited state is a sum of weighted phase 
factors in which the summation index m enters through 
the dimensionless offset ?7,„. Since according to Eq. ([1]) 
the offset is linear in m and the phase is quadratic in 
rjm the phase factors are quadratic in m. Such sums are 
called Gauss sums Il9l I20I . 



III. SINGLE INTERMEDIATE STATE 

In order to gain some insight into the dependence of 
the transition probability amplitude Cg on the dimension- 
less chirp parameter a we first concentrate on a single 



intermediate state |m). In this case the sum in Eq. pH)) 
reduces to a single term and the transition probability 
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|erfc(Cm(a))|' (29) 



is solely determined by the weight factor Wm- As a con- 
sequence the quadratic phase r^^a in Eq. (j26p has no 
chance to contribute. Nevertheless, under appropriate 
conditions 0, 0, H, [IBl we find interference due to the 
rather intricate properties of the complementary error 
function. In order to bring this subtle phenomenon to 
light we now analyze the asymptotic behavior of erfc. 



A. Path of integration induced by chirp 



According to Eq. the transition probability W^^ 
follows from the absolute value squared of the comple- 
mentary error function evaluated at the argument Cm(a)- 
The function erfc(Cm) is an integral of the analytic func- 
tion 



(30) 



along the path C(Cm)- For \zi\ < \zr\ the integrand 
exp(— 2;2) decreases rapidly for increasing \zr\. The de- 
caying domains for positive and negative real parts are 
connected via a saddle at the origin. Indeed, in the sec- 
tor \zr\ < \zi\ we deal with a rapidly growing integrand 
as \zi\ increases. 

As a function of the chirp parameter a the starting 
point Cm (fl) of the integration traverses a path in complex 
plane. When we square the representation 



Cm = Vmix + iy) 

and take the real and imaginary parts of 



y + 2i xy = i a ~ 1, 



we find the trajectory 



and the translation 



2xy = 2x \/l + 3p 



(31) 
(32) 

(33) 
(34) 



between x and a. 

For large positive or large negative values of x the tra- 
jectory, Eq. ([33|) . of Cm approaches the diagonals of the 
first or second quadrant, respectively. The straight path 
C of integration defined in Eq. (|2ip has a steepness de- 
termined by the coefficient 



■(l+a2)-l/4e5a'-<=tan(a) ^35) 



in front of t. Hence, for large positive values of a the 
straight path encloses an angle with respect to the pos- 
itive real axis which is slightly smaller than 7r/4. For 
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|erfc(Cm(a))| 



FIG. 2: Complementary error function in its dependence 
on the chirp parameter a (left bottom) as obtained by inte- 
grating along a straight path in complex space (top). The 
starting point C^m(a) of this path traverses as a function of 
a a trajectory in complex space which for large positive and 
large negative values of a approaches the diagonals of the first 
and the second quadrants, respectively. The path is always 
a straight line which is tangent to this trajectory and tends 
towards infinity as shown on the top. The gray curve below 
translates between a and Re(^m. 



1. Decaying regime 

We start our discussion by analyzing the behavior of 
erfc(Cm(a)) for large positive values of a. For this purpose 
we integrate the integral in the definition Eq. p4| by 
parts which yields 



erfc(Cm) = 



dz 



(36) 



Since the path Cj^ shown in Fig. [3] reaches the domain 
where decays the boundary term originating from 
the end point of the path C+, that is from infinity, van- 
ishes. As a consequence only the starting point of C_|_ 
contributes leading to 



erfc(Cm) 



1 



dz 



(37) 



We can now repeat this procedure of partial integration 
to obtain a power series of erfc in I /Cm- In lowest order 
we arrive at the asymptotic expression 



erfc (Cm) 



1 



(38) 



valid for 1 <^ a. 



2. Oscillatory regime 



a = the path is parallel to the positive real axis. For 
large negative values of a the angle of inclination of C is 
slightly larger than — 7r/4. 

In Fig. [2] we show on the top by a solid line the trajec- 
tory y{x) given by Eq. (|33p together with the direction 
of the path C indicated by arrows. Below we depict the 
translation function a{x) connecting the chirp parameter 
a with the real part of Cm- The picture on the bottom 
left shows the absolute value of the complementary error 
fimction erfc(Cm(a)) as a function of a. 



B. Asymptotic expressions 

The top of Fig. [T] and the bottom left of Fig. [2] demon- 
strate that the value a = of the chirp parameter sep- 
arates two distinct domains of the transition probabil- 
ity: (i) For positive values of a We decays as a increases, 
and (ii) for negative values of a the probability oscillates. 
These qualitatively different behaviors of We originate 
from the different slopes of the paths of integration as we 
now show. 



We now turn to the limit of large negative values of 
a. It is tempting to apply again the technique of integra- 
tion by parts. However, this idea is bound to fail since 
the partial integration creates inverse powers of z in the 
integrand. As shown in Fig. [2] for large negative values 
of a the path of integration gets very close to the origin 
of complex space. As a result the series cannot converge 
in this case. This heuristic argument clearly indicates 
that the behavior for negative values of a is substantially 
different from the one arising for < a. 

In order to evaluate erfc(Cm(a)) for large negative val- 
ues of a we first recall that is an analytic function. 
Therefore, the Cauchy theorem applied to the closed path 
Cc depicted in Fig. [3] yields the identity 

j)dze-^''^0. (39) 

Cc 

Since in infinity the integrand vanishes the contributions 
from the paths Coo and C_oo are zero and we find imme- 
diately 
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excitation process it is more convenient to use the repre- 
sentation Eq. (I7|) for the transition probability ampHtude 
Ce with the integral Im ■ The latter is the product of two 
integrals which reflects the sequential character of the ex- 
citation process: in the first step, the atom is promoted 
to the intermediate state |m), while the second step pro- 
vides the residual energy to reach the excited state. 



1. Times of excitation 



FIG. 3: Paths C+ and C_ of integration in the complex plane 
corresponding to a positive and a negative value of the chirp 
parameter a . We depict the domains in the complex plane 

_ 2 

where the integrand e ^ of the complementary error function 
increases by shaded areas. For < a the starting point C^m of 
the integration is located in the first quadrant and the straight 
path C+ is inclined with a steepness slightly smaller than the 
first diagonal. Thus C+ proceeds almost completely in the 
domain of the complex plane where the integrand exp(— a^) 
of erfc decays. For a < the path C_ of integration starts in 
the second quadrant and proceeds slightly above the second 

_ 2 

diagonal into the decaying region of e ^ . In order to evaluate 
the resulting integral we close the path as indicated by the 
dashed line consisting of the paths C_, Coo, Cr, C-oo and C+. 



The path C+ is the mirror image of the path C+ and 
avoids the origin. Therefore, we can use the technique 
of the preceding section to approximate the remaining 
integral. When we traverse C+ in the opposite direction 
we arrive at 



erfc(Cr„) = 2 + 



1 



(41) 



In this approximation the complementary error function 
is a sum of two complex- valued contributions: a constant 
term and a phase factor. As a result the interference 
between these two terms manifests itself in oscillations of 
the probability We as a function of a. 



In order to bring out this feature most clearly and to 
identify the times when the transitions occur we now sub- 
stitute the expression Eq. (fTT|) for the shape f{t) of the 
pulse into the definition Eq. ^ of X,„ which leads us to 
the integral 



dt' e 



-oo 

t' 



xe 



dt" e-"-*'"e-^"'(*"- 



(42) 



So far the calculation is exact. We now recognize that 
for large values of the chirp parameter a the definition 
Eq. (|15p of ar and provides us with the inequality 



ar 



1 

= ^ < 1. 
a 



(43) 



Consequently, the real-valued Gaussians in the integral 
Eq. p2|) arc slowly varying compared to the oscillatory 
terms resulting from the quadratic phases. The main 
contributions to the integral emerge from the oscillatory 
terms and in particular from the neighborhood of the 
times 



t'! 



2a, 



Vrn 
— C 



1 



1 



Aw 



(44) 



and 



IV. INTERFERENCE OF EXCITATION PATHS 

The preceding section has shown that the two paths 
C+ and Cr of integration in complex space correspond 
to two interfering contributions in the total probability 
amplitude Ce- We now show that they represent two 
different paths of excitation in the atom. Indeed, the 
term e"^™ /Cm results from a direct two-photon transition 
whereas the term 2 results from a sequential excitation. 



A. Sequential excitation 

So far our analysis was based on the integral I2 ■ How- 
ever, in order to identify the sequential character of the 



2a,; 



(45) 



when the phase factors are slowly varying. 

Therefore, the transition from the level \g) to the level 
|m) appears at the time t" ~ t" followed by the transition 
from |m) to |e) at the time t' = t'^. In order to ensure 
this time ordering enforced by the limits t" < t' of the 
integral Xm defined in Eq. (|42)) we need to have 



t'l < t' 



(46) 



This inequality is only satisfied provided is negative. 
From Eq. (|^ we recognize that it is the product 77™ a 
which determines the sign of t'^. Consequently we only 
find points of slowly varying phase for 77^0 < 0- In th^ 
case < rj^a no such points exist and the integral X„j is 
decreasing rapidly. 



2. Origin of quadratic phases 

Moreover, this analysis of the integral X,„ brings out 
most clearly that each transition is associated with a 
phase which is quadratic in the dimensionless offset rjm- 
The phases of both transitions are identical. The total 
acquired phase in the two-photon transition is twice of 
that of the individual one-photon transition. 

We can now evaluate the integrals in an approximate 
way by replacing the variables t' and t" in the real- valued 
Gaussians by t'^ and and factoring them out of the 
integral which yields 



A E/h 



















\9) , 




i 

1 



Im = &{-T]ma)f^ 



X exp 





du 






oo 






exp 


US.) 


ai 2ai 


V 2a J 



(47) 



Here we have also extended the integration over t" to 
-l-oo. 

With the help of the integral relation [2^ 



du f 



(48) 



we find 



l-m ~ Q{-il7na)fo - — exp 
I a,. 



52 
ai 2a,; 



a^ 



exp 



2a,: 



(49) 

When we now take the limit of 1 ^ |a| in the definitions 
Eq. ((T2|) and Eq. ([15]) of /o and a^ which yields f§ =i/a 
and ai — Alu'^ /{2a) we obtain together with ar/ai — 1/a 
the final approximate expression 



27T 



e{-7]^a) exp {-r|l^) exp {i tjI^o) 



(50) 



We substitute this formula into Eq. ([7]) for the probability 
amplitude Ce and arrive at Eq. ([26]) with the approximate 
weight factor 



= -c?m0(-'7ma) exp (-7?^) 



(51) 



This expression also results from the exact formula 
Eq. (H?]) with the help of the relation VI - ia = (-1)^/^0 
together with the primitive asymptotic expansion 

erfc(C,„)«2e(-77,„a) (52) 

of the complementary error function following from 
Eqs. (pll and (jiT]). 

3. Interference m time-frequency phase space 

This analysis of the transition probability amplitude in 
second order perturbation theory translates itself into an 



FIG. 4: Origin of Gauss sum from phase space interpreta- 
tion of two-photon excitation induced by a chirped pulse. In 
the phase space spanned by frequency v and time t we in- 
dicate the ground state \g), the intermediate level |m) and 
the excited state |e) of frequencies 0, uq + 5m and 2u}o, re- 
spectively, by thick horizontal lines. We also denote the one- 
photon resonance of energy luq and the virtual level of fre- 
quency ujo — Sm by dotted horizontal hues. The instantaneous 
frequency uj{t) — iJo + 2a.i t is depicted for one negative value 
of the chirp by the tilted thin line. A stationary phase analy- 
sis of the double integral Eq. (|8} determining the two-photon 
transition probability amplitude brings out most clearly that 
the sequential excitation appears at the times = 5m/(2ai) 
and t'g. Indeed, at the instantaneous frequency is resonant 
with the transition from \g) — > |m) whereas at t!^ it is resonant 
with the transition from |m) to |e) as indicated by the dotted 
and the thick arrow. In this phase space representation the 
times when the transitions occur are determined by the cross- 
ings of the tilted thin line with the thick and the lower dotted 
horizontal lines corresponding to the frequencies loq -\- (5m and 
— 5m- The enclosed shaded area Ai + A2 corresponds to 
the phase acquired during the transition. Each triangle has 
the area 5mta/2 = 5^/(40^). As a consequence, the phase 
associated with this transition is quadratic in the offset Sm- 
In the case of a multi-state ladder with a harmonic manifold 
we have the interference of all excitation paths, that is the 
weighted sum of phase factors with quadratic phases giving 
rise to a Gauss sum. 



elementary geometrical representation in a phase space 
spanned by frequency v and time t. In Fig.|4]we consider 
a single sequential excitation path in this space. We in- 
dicate the frequencies 0, t^o + Sm and 2luo of the ground 
state, the m-th intermediate state and the excited state 
by thick horizontal lines. We also denote the frequency 
luq of the two-photon resonance and the virtual level of 
frequency loq — Sm by dotted horizontal lines. The in- 
stantaneous frequency Lu{t) = ujq + ait is depicted for 
one negative value of the chirp by the tilted thin line. 

The transitions from \g) to |m) and from |m) to |e) 
occur at times t" and t'^, respectively. These times are 
determined geometrically by the crossings of the repre- 
sentations of the instantaneous frequency uj(t) and the 
two frequencies loq -\- Sm and LjQ — Sm-, that is by the cross- 
ings of the tilted thin line with the thick and the lower 
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dotted horizontal lines. Each crossing contributes to the 
total transition probability amplitude with a phase de- 
termined by the enclosed triangular shaded area Ai + A2 
corresponds to the phase acquired during the transition. 
Each triangle has the area t"(5,„/2 = 6f„/{2ai). As a 
consequence, the phase associated with this transition is 
quadratic in the offset 5,„. 



B. Direct excitation 

We conclude by briefly discussing the direct two- 
photon excitation. In the analysis of the time ordered 
integrals Eq. ((42|) we have implicitly assumed that the 
moments t'^ and <" of the transitions are appropriately 
separated in time. However, since we deal with a dou- 
ble integral the instant t' = t" can make a significant 
contribution. This fact stands out most clearly in the 
re- formulation, Eq. of the double integral where the 
line t' = t" translates into the condition t = 0. When 
we recall the definition Eq. (|2T|l of the path C(C) of in- 
tegration the moment i = translates into the starting 
point C,rn of the path. As a consequence the main con- 
tribution to the integral I2 arises from the lower limit. 
The partial integration that we have used in Sec. lIIIBl to 
derive the asymptotic expansion of the complementary 
error function took advantage of this fact. 



V. MANIFOLD OF INTERMEDIATE STATES 

So far we have concentrated on a single intermediate 
state. We now briefly address the problem of a manifold 
of intermediate states. In the case of a multi-state ladder 
with a harmonic manifold we have the interference of 
all excitation paths, that is the weighted sum of phase 
factors with quadratic phases giving rise to a Gauss sum. 
In the case of a positive chirp parameter a these paths 
consist solely of direct two-photon excitations. In the 
limit of negative values of a each interfering path through 
the m-th level consists of an interference of a sequential 
and a direct excitation. However, the direct excitation 
is less pronounced especially for appropriately negative 
values of the chirp. In this case the probability amplitude 
is predominantly determined by the interference of the 
quadratic phases arising from the sequential excitations. 



dispersion a to verify the quadratic nature of the phase 



5m 

Alu 



a = -N£_ + 27r ( m + — ) e (53) 



N 



in the probability amplitude Cg. Here we have introduced 
the quantity 



N = 

and the dimensionless chirp 
6A a 



25 
A 



5A 



TT AlU^ 



(54) 



(55) 



Moreover, in terms of N the offset 6m defined by Eq. ([T]) 
reads 



A ( m+ ^ 



A{ m + — 



(56) 



Hence, the final expression for the excited state proba- 
bility amplitude takes the canonical form 



M 

exp (^i ^A^C) X! ^™ *^^P 



( m + — ) ^ 



(57) 

of a Gauss sum. Here we have numbered the intermediate 
levels from —M' to M. 

When we express dm with the help of Eq. ([56]) the 
weight factors 



Wm = — 1^ erfc (Cm) exp 

as well as the argument 

TO -I- A^/2 . 



m + N/2 
Am 



Cn 



Am 



(58) 



(59) 



of the complementary error function contain the width 
Am = Aw/A. 

The transition probability amplitude consists of a sum 
of weighted phase factors with phases that depend lin- 
early and quadratically on the summation index to. 
Closely related sums appear in the context of fractional 
revivals [H [Hi, the Talbot effect [111 and curlicues [2i,[23| 
and Josephson iunctions(28|. 



A. Gauss sum made explicit 

In order to bring out most clearly the connection to 
the Gauss sums, we make the dependence of the weight 
factors and the phase on to more explicit. For this pur- 
pose we use the definitions, Eq. ^ and Eq. for the 
offset Sm of the harmonic manifold and the second order 



B. Asymptotic expansion of weight factors 

With the help of the asymptotic expansion Eq. ((5T|) 
the approximate real- valued weight factors 



-dm. 6 



N 
1 



exp 



7V/2 



Ar 



(60) 
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are given by half a Gaussian centered around m w —N/2. 
The sign of the chirp determines which half. Indeed, for 
negative chirps that is a < we find only contributions 
for TO > —N/2. In contrast for positive chirps, that is 
for < a only levels with m < —N/2 contribute. More- 
over, by controlling the width Am we can ensure that 
all quantum paths enter the sum with approximately the 
same weight. 

In Fig. [5] we compare and contrast the exact (black 
dots) and the approximate (gray dots) expressions, 
Eqs. ([SS]) and ([60]), for the ratio \wm/dm\ of the weight 
factors Wjn and the matrix element dm- To be specific, 
we have used iV = 15 and a negative chirp. The figure 
clearly shows that in this case the weight factors are non- 
vanishing only on the right half of the Gaussian, that is 
for —N/2 = —7.5 < m. Moreover, there is an excellent 
agreement between the exact result and the approxima- 
tion. 

We conclude by noting that there are small oscilla- 
tions in the exact weight factor distribution which are 
not contained in the approximate expression. They are 
a remnant of the interference between the direct and the 
sequential path as discussed in Sec. IIVI 
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FIG. 5: Comparison between the exact (black dots) and 
the approximate (gray dots) expression for the scaled weight 
factors |wm/dm| given by Eqs. (|58|l and (|60|l . respectively. 
Here we consider the case A'' = 15 with Am — Auj/ A = 50.8 
for the negative dimensionless chirp a — —10824 < 0. 



VI. CONCLUSION AND OUTLOOK 

In the present paper we have analyzed a two-photon 
transition in a ladder system driven by a chirped laser 
pulse. We have shown that the resulting excitation prob- 
ability amplitude is of the form of a Gauss sum. However, 
one task remains: Find a quantum system which displays 
such an equidistant ladder system. 

Before we address this task it is helpful to recall the 
origin of the Gauss sum - the interference of paths with 
quadratic phases. Fig. |3] brings out most clearly that 



these phases arise already from the first step of the exci- 
tation, that is from transition into the equidistant man- 
ifold. Therefore, a one-photon transition from a ground 
state to an equidistant manifold of excited states driven 
by a chirped pulse also leads to a Gauss sum. In this 
simplified system the fluorescence signal is acquired from 
the states in the harmonic manifold. Interference from 
the individual levels occurs during spontaneous emission 
of the states. 

At least five different quantum systems offer them- 
selves as candidates for such one-photon excitations: (i) 
Rydberg atoms in electric or magnetic fields, (ii) vibra- 
tional energy levels in a molecule, (m) atoms in traps, 
(iv) laser driven one-photon transitions, and (v) quan- 
tum dots. However, these systems have to be subjected 
to additional requirements: (i) the Rabi-frequencies for 
the interfering excitation paths should ideally be of the 
same order of magnitude in order to guarantee that they 
contribute in a democratic way. that is with the same 
weight, (ii) since anharmonicities in the target manifold 
change the character of the resulting probability ampli- 
tude, only perfectly equidistant states result in a Gauss 
sum. We now briefly discuss these quantum systems and 
address these requirements 

Rydberg atoms [29| are highly sensitive to electric or 
magnetic fields. In an electric field the energy levels of 
a Rydberg atom split into a manifold of statesjs^. The 
so-called Stark maps which show the dependence of the 
splitted spectral lines as a function of the electric field 
have been measured for many alkali atoms. Whereas for 
a modest field strength the separation between neighbor- 
ing splitted energy levels is constant, for larger values 
of the field also higher order corrections to the energy 
have to be taken into account. A closer analysis of the 
Stark map reveals that the dipole moments vary for each 
member of the manifold and the individual paths do not 
contribute equally. Due to the Zeeman effect a homoge- 
neous magnetic field also creates a manifold of equidis- 
tant states \m). However, in this case the problem are 
the selection rules which eliminate most of the transi- 
tions. A change of the orientation of the magnetic field 
might offer a possibility to overcome this problem. 

Diatomic molecules may also serve as a candidate 
system for identi fying a harmonic manifold of (almost) 
equidistant states [3l|. The diagonalization of the ro- 
vibrational Hamiltonian yields that the inter-nuclear po- 
tential in a diatomic molecule can be indeed approxi- 
mated as a harmonic oscillator with equidistant spec- 
trum. However, anharmonicities have again to be taken 
into account and would destroy the regularity of the 
Gauss sum. 

Another promising realization of an equidistant spec- 
trum relies on cold trapped atoms or ions [32,, The 
central idea of this approach is that the trapping in mag- 
netic microtraps relies on the magnetic dipole of the 
atom which depends on its internal state. Therefore, the 
center-of-mass motion of an atom depends on its internal 
state. Our suggestion for the realization of an equidistant 
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spectrum is reminiscent of the quantum tweezer |34l |35|- 
Whereas the atom initiaUy in the ground state is weakly 
bound to a shallow trap, the atom in its excited state feels 
a steep harmonic potential. The vibrational states in the 
excited electronic state constitute the desired harmonic 
manifold. 

Yet another method to obtain an equidistant ladder 
system is to apply two electromagnetic fields to a two- 
level atom which has a permanent dipole moment in the 
excited state. A strong cw-modulation field generates 
an equidistant Floquet ladder in the excited state. The 
second field, that is the weak chirped pulse, leads to the 
excitation probability amplitude in the form of a Gauss 
sum. 

Quantum dotsjsHl may be the perfect choice in the 
end. Indeed, they allow for designing the discrete en- 
ergy spectrum|37| of a trapped electron and can thus be 
viewed as an artificial atom. Quantum dots have already 
been discussed as a candidate system for implementing 
quantum logic [ssj. 

In conclusion, we have shown that a two-photon tran- 
sition in a ladder system driven by a chirped laser pulse 
leads to an excitation probability amplitude in the form 



of a Gauss sum. We have discussed several physical sys- 
tems in the context of realizing such an excitation scheme 
and have also addressed the experimental challenges. We 
emphasize that this example also establishes an interest- 
ing connection between number theory and wavepacket 
dynamics. 
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